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Abstract 

The most useful measure of a bipartite entanglement is the von Neumann entropy of either of 
the reduced density matrices. For a particular class of continuous-variable states, the Gaussian 
states, the entropy of entanglement can be expressed rather elegantly in terms of the symplectic 
eigenvalues, elements that characterize a Gaussian state and depend on the correlations of the 
canonical variables. We give a pedagogical step-by-step derivation of this result and provide some 
insights that can be useful in practical calculations. 
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I. INTRODUCTION 

Entanglement is recognized as a key resource in quantum computation A num- 
ber of protocols utilize entanglement for performing tasks that would be very hard, if not 
impossible, in classical information processing j^J. Different quantitative measures cap- 
ture different aspects of entanglement 3j. For pure states there is a preferred way to 
quantify the entanglement between a subsystem A and its complement B. This is the 



degree of entanglement or entropy of entanglement [4], defined as the von Neumann en- 
tropy of either of the reduced density operators pa or p B of the state p = pab- Explicitly, 
S = — tr (p^logp^) = ~ tr (p^logps). While studies of entanglement originally focused on 
finite-dimensional quantum systems, continuous- variable systems are becoming increasingly 



important both practically and theoretically js|. For instance, a special class of continu- 
ous variable states, the Gaussian states, plays an important role in quantum optics and 
quantum information processing and can be described conveniently by an easy algebraic for- 
malism. Furthermore, entanglement calculations in infinite-dimensional systems are much 
more demanding than their finite-dimensions counterparts. Unlike the general case, the von 
Neumann entropy of Gaussian states has a simple expression in terms of a finite number 



of the symp 
Specifically 



ectic eigenvalues {<7j} of the correlation matrix of the quadrature operators. 
6], 

n sub 



S(p) = i>2 (*i + !og 2 (*i + ^ - (*i - ^ log 2 (ffi - 



where the index % runs over the modes of one of the two subsystems under examination. In 
this paper we give a simple but precise derivation of this formula, showing the link between 
covariance matrix of a Gaussian state, symplectic eigenvalues and entropy. First we present 
a short introduction to Gaussian states, defining a covariance matrix and its symplectic 
eigenvalues, and then we give a derivation of (0Q). We conclude with some observations on 
entropy and temperature. 

II. GAUSSIAN STATES 

Gaussian states are the basic ingredient of continuous-variable quantum computation 
[7, 8]. They are well- understood, easy to produce in a laboratory and can be described using 
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a simple matrix formalism. Each mode of a quantized electromagnetic field is equivalent 
to a quantum harmonic oscillator and the quadrature operators q and p are related to the 
mode creation and annihilation operators by jo]: 



(a^ + a), p = \/hm(jj—=(a) — a) . (2) 

Consider a bosonic system having 2n canonical degrees of freedom (such as n light modes), 
and introduce a 2n-dimensional column vector r = ...,q n ,pi, ...,p n ) T of quadrature op- 
erators. The canonical commutation relations 

[ft, Pj] = , [ft, Qj] = , [Pi, Pj] = , (3) 
can be rewritten in matrix form using the components of r as 



htt itj = -t[r h rj] , = [ n \ . (4) 




The skew-symmetric matrix Q is called symplectic (metric) matrix 10|, Notice that 
Q T = Q- 1 = -Q. 

The ground state of a quantum harmonic oscillator is a Gaussian state. It is clear intuitively 
from the Gaussian shape of its wave function and will be given a precise meaning shortly. 
Similarly, the ground state of a system of n harmonic oscillators, whose Hilbert space is 
given by the tensor product of the individual Hilbert spaces % = (S)"=i ^ s a Gaussian 
state. In fact, the ground state of any system that is described by a Hamiltonian quadratic 
in the canonical operators, 

1 n n 

i j 



defined by a real and positive-semidefinite crossing matrix H, is Gaussian [12|. Mathemati- 
cally this means that the state characteristic function, which we are going to define now, is 
a Gaussian in the phase space 

Any state p of n modes can be characterized by the following construction on the phase 
space of the system. We introduce the Weyl operator 



n n 



W v = exp{-t V T Qr} = (g) A(«i) = (g) e*^-^ , (6) 

i=l i=l 

which in quantum optics is taken as a phase space displacement operator. Here the vector 
r] = (ai, a n , bi, b n ) T e M 2n defines the displacement that can be represented in the 



complex form as cc, = ^(oj + ibi) G C. The action of the i-th mode displacement operator 



Di(a.j) on the i-th mode ground state results in the coherent state \ai) = Di(oti)\0) 



The characteristic function of the state p is defined as the expectation value of the Weyl 



operator (for a textbook reference see for example 151]). 

X P {rf) = tr[ P W ri ] . (7) 

It is equivalent to the Wigner distribution, which is defined in terms of phase-space variables 

as 

W{q, V ) = ^f d n q'{q-q'\ P \q + q')e 2iq ' p , (8) 

and is a phase space representation of the density matrix p. The Wigner function is usually 
expressed by the symplectic Fourier transform of the characteristic function, 

W(X) = ^ J d^ V e^ x xM . (9) 

with X = (gi, ...,q n ,pi, ...,p n ) T . The density operator of the quantum state can be written 
in terms of its characteristic function by means of a Fourier-Weyl relation 



P = WrJ d VXei-vW^ (10) 

where W(rj) is the Weyl operator, and therefore the state is uniquely determined by its 
characteristic function \p- Finally, a state p of n modes is said to be Gaussian whenever the 



characteristic function has a Gaussian shape in phase space 16[], which means that it can 
be written as 

XM = Xp (0)e~^™ T ^ DT ^ . (11) 

The first two statistical moments are captured by the vector D = tr[prj] = (rj) of expectation 
values of the quadrature operators and the 2n x 2n real symmetric matrix T that carries the 
information about the variances 

T i)j = -RBtt\p(r i -{r i ))(r j -(r j ))] . (12) 

All higher-order statistical moments of a Gaussian state can be expressed from D and T. 
The correlation matrix V is called covariance matrix (or sometimes noise matrix) and it plays 



a central role in the fo 



change entanglement [3] and hence, since displacements are single modes local translation in 



lowing discussion of entropy. Local unitary transformations do not 



phase space, they leave the entanglement properties of the state unaffected [171 ] . Therefore 
the elements of D contribute nothing to the entanglement and they can all be made zero. 
Accordingly, we can rewrite the covariance matrix as 

Tij = Retr [pr^r,] . (13) 

However, the matrix T cannot be arbitrary. In order for a real symmetric matrix to be 
the covariance matrix of some physical state we need to fix a constraint: The canonical 
commutation relations require the positive definiteness of 

r + izO^O, (14) 

which is just another way to rewrite the Heisenberg uncertainty relations fl8 |. 
Gaussian states are important because of the existence of a class of operators corresponding 
to common laboratory procedures that preserve the Gaussian properties of the states on 
which they act. Thus we define a Gaussian unitary operation as a transformation that maps 
a Gaussian state onto a Gaussian state. 

There exists a symplectic representation of the Gaussian unitary group lfjj. To each Gaus- 
sian transformation U we can associate a unique symplectic transformation S G Sp(2n, R) 
(for more details about the real symplectic group and its properties see 10]). The group 
element 5" describes a linear transformation of the quadrature operators expressed by 

p< = U(S) p U\S) — > r' = Sr = UiS^rUiS) . (15) 

These transformations preserve the commutation relations. Therefore, the action of any 
S G Sp(2n, R) on the matrix Q is given by 

in = [rl, r'j] = S[ri, r y ]£ T — ► Q = SQS T . (16) 

Furthermore, if S is a symplectic transformation then it also satisfies S T = QS^ 1 ^ 1 , S^ 1 = 
nS 7 ^ 1 G Sp(2n, R). Under the action of a symplectic transformation the covariance matrix 



transforms as 



111 

T' = cov (Sr) = S cov(r) S T = STS T . (17) 



A. Symplectic eigenvalues 



The next step towards understanding of the entanglement entropy of Gaussian states is 
to introduce the concept of symplectic eigenvalues. In fact, equation ffiTj) tells us that the 

5 



covariance matrix does not transform by a similarity transformation under the action of 
the symplectic group. We are therefore interested in finding an alternative form of T such 
that its eigenvalues are invariant under a symplectic transformation and can thus uniquely 



characterize the state. To do so we use Williamson's theorem [19[. It states that any real 
symmetric positive-definite 2n x 2n matrix, such as the covariance matrix T, can always be 
made diagonal by means of a suitable symplectic transformation S w G Sp(2n, M), 

S W TS W = r^, , (18) 

where now 

T w = diag(cxi, a 2 , cx n , cri, ct 2 , a n ) (19) 

and all the <Tj are real. Hence, after the transformation S w we have, for the transformed 
canonical operators r', that Re(f^f^) = 5ij<Ji. The matrix T w is called the Williamson 
normal form of the matrix T. It is important to realize that in general the {cr^} are not the 
eigenvalues of T or of any T$ = S T TS determined by a transformation S different from S w . 
Now define a new matrix M such that T = —MVt and thus 

nr 1 = -Afiwr 1 m = m . (20) 

Using (Tl6l) it is easy to see that S T Q = QS" 1 and therefore the action of a symplectic 
transformation S on the matrix M results in a similarity transformation of M that preserves 
its eigenvalues, 

m' = t'q = srs T rt = sms- 1 = sms- 1 . (21) 

Hence every matrix M determined by varying S over the group Sp(2n, R) shares the same 



spectrum [161 ] . In particular, if we take the matrix T to its Williamson form and transform M 
accordingly, the eigenvalues of the matrix M' = S w TS^,fl = T W Q will be equal to {±ia{} [20 ] . 



The n absolute values {cr^} of the elements of the spectrum correspond to the n dinstinct 
eigenvalues of 1TV We say that these are the symplectic eigenvalues of the (covariance) 
matrix V and call the set {cij} the symplectic spectrum. This spectrum characterizes the 
Gaussian state. By construction it is invariant under any Gaussian transformation. 
Once the covariance matrix V of the state is given or calculated, the set of symplectic 
eigenvalues can be directly obtained from the spectrum of the matrix M = TQ. We will 
see in the next section how the symplectic eigenvalues contain the total information about 
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entanglement properties of a Gaussian state. This is the reason why this algebraic description 
of Gaussian states is efficient. It allows one to quantify entanglement simply from the 
symplectic eigenvalues of the matrix of correlations of the quadrature operators, which are 
in general much easier to calculate than the eigenvalues of the density matrix. 



III. DERIVATION OF THE FORMULA 



In this section we give a step-by-step derivation of the formula ([T]) for the von Neumann 
entropy of a Gaussian state, showing explicitly the connection between symplectic eigenval- 
ues of the covariance matrix and entropy. 

Consider a general Gaussian state p = ^2jPj\<f>j){(j)j\ corresponding to n modes with 
a covariance matrix T. We can always find a symplectic transformation S such that 
p' = U(S) pU^S) — > f' = Sr, which takes T to the normal form T' = STS T . 

Note that in general this is not the Williamson form. Hence, we can rewrite p' as 

p' = p' l ®p' 2 ...®p' n , (22) 



where each p\ is the density matrix of a single thermal oscillator 23| . Observe that after the 
transformation, the transformed oscillators are now uncoupled non-local thermal oscillators. 
A harmonic oscillator in thermal equilibrium at temperature T is described by a canonical 



ensemble Jq]. In the number basis {| </?„)}, its density matrix is written as [24] 



P = ^2Pn\tPn){tPn\ , (23) 
n 

where the p n 's are the probabilities associated to each state and correspond to p n = 
Z v e~ En l kBT ^ with E n energies of the n-th state of the Hamiltonian H\(p n ) = E n \(p n ), and 
where the partition function Z = Tr ^e~^/ kBT ^j is a normalization constant. Hence, the 
density matrix of each uncoupled mode is equal to 

n 

= zr 1 e- 6 V kBT * . (25) 
Using the transformed creation and annihilation operators 
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whose transformation under S follows from that for quadrature operators, 



£ = (&!,... ,a n ,a\,...,al^ — > £' = S£ 



(27) 



we can rewrite the Hamiltonian of each oscillator as H[ = huu'^a/a^ + |) and the partition 



function Z$ as 

oo 



t '^+|)K/*BT 4 



(28) 



n' i= n' 1= 

where v! i is the eigenvalue of the number operator at'c^|<p n /) = n^|</? n <) and = hw'JksTi is 
a cumulative parameter that depends on the transformed frequency u[. We can rewrite the 
last bit using the properties of the geometric series 

e -ft/2 



-n'A 



2v 



n.'=0 



1-e- 



(29) 



and finally we have 

= (1 - e -ft) e-al'^A . (30) 

To simplify this expression we can rewrite the density matrix in terms of = (n'J, the 
mean occupation number of the transformed modes: 



< = <alX> = Tr(P^X) = J> n ,| (1 - e" ft ) e-^ai'af^ 



n'=0 



1 + n' 



n'=0 



e ft _ i 

As a result we obtain the following expression for the density matrix 

1 



(31) 



/ 

Pi 



1 + n' V 1 + n 



— i \ a a- 



n 



(32) 



Calculation of the von Neumann entropy 5*(p^) = — Tr (p' { logp'J gives: 

at'*/ 



n'=0 



log 



n =0 



log(l + n-) + n- log 



IVn<) 



'/I: 



1 + n' 



IVn<) 



1 + n' 



M 1+ ,DE(r^)"' + '-(r^)E<T^)"' 

n =0 V J/ V l/ n =0 V l/ 



logn'i - (1 + n-) log(l + n-) 



(33) 
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where we used 

oo 

E(lfH;j ' = (i + ^ E^TTiiJ "=^ 1 + ^ (34) 

It immediately follows that for a single oscillator thermal state, the von Neumann entropy 
is expressed, in terms of mean occupation number, as 

S{pft = (1 + n[) log(l + nj) - n[ log 7^ , (35) 



a well-known result in statistical physics (24]. 



The connection of the entropy with the symplectic eigenvalue a\ of the state follows straight- 
forwardly. After the symplectic transformation on the canonical variables that takes T to 
its normal form, the covariance matrix of the reduced state p\ looks like 

r; = ( &) ) . (36) 

Using the ladder operators f[2T)j) it is easy to show that 

m = Trm') =z-> E^i^y^' + + ^e-m+m/^^ = (37) 

1 h cotlA (38) 



2m[w; m'jOu'i - 1 2m^ 2 
and, proceeding in the same way, that 

(p?') = — - + hrdw'j-s = — — coth — . (39) 

Our previous discussion allows us to exhibit a relationship between the symplectic eigenvalue 
a\ of the system and the mean occupation number. We know that the spectrum of I^fi 
corresponds to {±2<r-}, and therefore from ( I3TT) we have that: 



■ fte ft + 1 , 1. 



eigenvalues {r'^} = ±iy/{%' ){pf) = ±%^— 7 = + o) (40) 

and after fixing the units such that h — 1, we find the equivalence 

^ = ^ + 5— >*J = ^-^ (41) 

Entropy is an additive quantity, therefore the total entropy of a state p' which is the direct 
tensor product of n states is just the sum of the entropies of each state. It is easy to rewrite 



the formula (135]) as a function of the symplectic eigenvalues, 

n 

s( P [) = £ 

Remind that the symplectic eigenvalues are the invariants of the correlation matrix. Since 
the entropy of a Gaussian state is solely a function of the symplectic eigenvalues, entropy 
itself is invariant under a symplectic transformation. This means that S(p) = S(p'). Hence 
we can drop the prime in the formula ( l4*2l) and find (CD), which concludes the derivation. 
We also notice from equations (|3~T!) and (T4T1) that the thermal parameter of each oscillator 
depends on the correspondent symplectic eigenvalue: 

IV. CONCLUSIONS 



< + I) log 2 U + \ 



< ~ I) log 2 U - \ 



(42) 



We conclude with an example and a couple of considerations about entanglement and its 
thermal properties. 



Systems of coupled harmonic oscillators, like harmonic chains [2l|, 122J, are relatively well- 
understood and easy to describe mathematically. There exists an extensive literature about 
quantifying the bipartite entanglement of such systems. Among other reasons, this is the 
case because some of these entanglement measures may follow area laws jfj]. States of coupled 
harmonic oscillators are also interesting because they can exhibit Gaussian properties. In 
particular, when the interaction between the modes is quadratic, the ground state of the 
system is Gaussian. 

For a pure state the entanglement entropy and the von Neumann entropy coincide 
Thus we can use the formalism introduced earlier to describe the entanglement properties 
of an modes pure Gaussian state p with covariance matrix V . We have seen that through 
a proper global symplectic transformation S, the state p can be decomposed as the tensor 
product of n single thermal oscillator ( |22l) . Now the new non-local states p\ have to be pure 
as well for the properties of decomposition of a pure state. Hence the decoupled oscillator 
are all in their ground state and we assign to each of them a virtual temperature Tj = 0. 
From equation (T4"T]) it follows that the symplectic eigenvalues are all equal to <7j = \ and 
making use of ([!]) we correctly find that the total entropy of any pure state p is zero. 
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Divide now the n modes into two sets A = {A\, A a } and B = {Bi, Bb} such that 
a + b = n and hand them to Alice and Bob (this very last step is not strictly necessary, but 
it is always fun). In order to calculate the entanglement entropy S(pa) = S(pb) we need to 
obtain the symplectic eigenvalues that belong to one of the two partitions. To start, assume 
that Alice wants to study her part of the system. In general, after the division, the reduced 
density matrix pa = tr bP corresponds to a mixed state. Therefore a local decomposition by 
means of a local symplectic transformation Sa, 

PA = PlJV ® Pr,T r ® Pr+l,G r+1 -- ® Pa,G a ( 44 ) 

contains both thermal and ground state oscillators. 

Alice's symplectic spectrum {ai^.^A is obtained from the spectrum of Ta^a where Ta = 
tr^T is the reduced covariance matrix of the set A. Then Alice can use her symplectic 
eigenvalues into the von Neumann entropy formula for Gaussian states ([T|) and quantify the 
bipartite entanglement. Suppose that Alice obtained s symplectic eigenvalues satisfying 

*W > \ , (45) 

and a — s symplectic eigenvalues (J s -\-i,...,a — \- This means that the local decomposition ( 144)) 
of Alice's set of modes can be rewritten as 

P'a = p'l, Tl - ® P' S ,T S ® P's+1,G S+1 - ® Pa,G a • ( 46 ) 

Physically this corresponds to having s transformed thermal oscillators with thermal pa- 
rameter fa given by equation (j43p and a — s ground states oscillators. Notice that only the 
thermal oscillators contribute to the bipartite entanglement entropy. 



Botero and Reznik made this construction rigorous. They proved in [23( that, after identi- 
fying the two sets A = {Ai, A a } and B = {B±, 5^}, it is always possible to write the 
Gaussian pure state \(f>) = \4>)ab as 



I A.B 



v) AlBl ® ... ® \<p s ) AsBs ® \0) As+i a ® |0) 5 , (47) 



where A = {Ai, ...,A a } and B = {Bi, Bb} are the transformed modes resulting from 
the application of local symplectic transformations on the set A and B, and s is equal to 
the number of symplectic eigenvalues associated to Alice's reduced covariance matrix Ta or 
to Bob's one T B . It means that the state \4>)a,b can be rewritten as the direct sum of s 
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two-mode squeezed states, where each mode belongs to a different partition of the system, 
and n — 2s oscillator ground states. 

Each two-mode squeezed state Ify)^ § is given by the expression, 

= -^£ e "* n/ Vh>>i% , (48) 

V J n 

and the squeezing parameter f3j corresponds to the thermal parameter of the j-th thermal 
oscillator of the local normal-modes decomposition of pa (or 

We hope this dicussion will help to clarify the necessary steps to derive the entropy expression 
for Gaussian states and offer at the same time an extensive literature where the reader can 
find more details on the subject. 



A. Two coupled harmonic oscillators 



We want to present an easy case in order to show some explicit calculations. Imagine to 
have a system composed of two quantum harmonic oscillators with mass m and frequency 
u coupled in position and described by the following Hamiltonian: 



H 



mu 



(Pi + Pi) + + + A (<?i - hY 



(49) 



2rn JZ/ 2 

The ground state of this system is Gaussian. We want to calculate the bipartite entanglement 
between the two oscillators for the ground state. The global symplectic transformation S, 
described by the matrix 

(l -1 \ 
110 
1-1 
\0 1 1 ) 

gives two uncoupled oscillators with the new frequencies 



V2 



(50) 



1 — U! , 

mui z 



1 + 



4A 



(51) 



The corresponding covariance matrix for the normal modes, see ( |36 

\ 



r 



1 

7! 



( ^ 

2mu) 




o 



2mwa 







u 2 









(52) 



12 



can be transformed back using ([IT]) into the covariance matrix T of the system 

\ 



-lp// cT\-l 



r = s-t\s 



1 

7! 



/ l+Q 1-a 

Amau) Amau) 

1—a 1+q 

Amau) Amau) 






\m{l + a)uj \m{— 1 + a)u 



(53) 





\ im(-l + a)u im(l + «)cj ) 

It is straightforward to trace out the complementary degrees of freedom and obtain the 
reduced covariance matrix for the first (or equivalently, the second) oscillator 



i\ = r 2 



l+q 

Amau) 







lm(l + a)u 



(54) 



The last step is to calculate the only symplectic eigenvalue of the reduced covariance matrix 
from the spectrum of fl2x2^i,2, which is given by 

1 + a 



4v^ ' 



(55) 



and finally it can be used to quantify the bipartite entanglement using the bipartite entan- 
glement formula (CQ). 
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